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Abstract: We show how and why the short distance ("hard") interaction, which is 
calculated in perturbative QCD, provides a mass cutoff in Gribov's formula for photon- 
proton collisions. This enables us to find a new and more restrictive unitarity bound for 
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this process, cr(7*p) < C{ln-)2. We develop a simple model that consists of "soft" and 
"hard" contributions, which yields a qualitative description of the published experimental 
data over a wide range of photon virtualities (Q^) and energies {W). This model provides a 
quantitative way of evaluating the relative rate of the short and long distance contributions, 
in the different kinematic regions. The main results of the analysis are (i) that even at 
= and high energies the short distance contribution is not small, and it provides a 
possible explanation for the experimental rise of the photoproduction cross section; and (ii) 
at large values of Q^, the long distance processes still contribute to the total cross section. 
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1 Introduction. 



The total cross section of a hadron-hadron interaction is bound by the Froissart - Martin 
hmit 

(Ttot < C In^ - , (1) 
So 

where C = depends on the mass of the hghtest particle exchanged in the crossed channel. 
The bound is a consequence of s-channel unitarity, analyticity and crossing symmetry. In 
spite of the ambiguity in the determination of C, we suggest that at the presently available 
hadron accelerator energies one should look for phenomena associated with s - channel 
unitarity, rather than the absolute bound. Indeed, a careful study shows that the scale of 
saturation of s-channel unitarity in elastic pp reactions is above the Tevatron energy, while 
the saturation scale for diffractive channels is considerably lower, appearing at ISR energies. 
This qualitative theoretical study is strongly supported by the experimental observation that, 
whereas aei/cTtot grows all through the ISR - Tevatron range, the ratio (Jdifjlotot decreases 
with energy [Q. 

The study of unitarity and the Froissart - Martin bounds in DIS are more complicated 
than for the hadron - hadron case. These complications originate from ambiguities in the 
implementation of the unitarity constraints due to electromagnetic photon coupling and the 
absence of a proper elastic channel, as well as the introduction of the mass of the virtual 
photon as an additional kinematic variable. From a phenomenological point of view, we 
have to take into account the "soft" and "hard" , or alternatively, the long distance and short 
distance phenomena, as contributors to the total 7*p cross section, or the target structure 
function. This differs from the usual picture of the hadron - hadron collision, where the 
incoming time like particle masses are fixed, and the total cross section is determined by the 
"soft" ( long distance ) Pomeron. 

A remarkable simplification of the DIS analysis has been suggested by Gribov |^ in the 
context of DIS on a nuclear target. Gribov's main observation was that at high energies, the 
7* fluctuates into a hadronic system ( i.e. qq to the lowest order ) with a coherence length, 

= — , which is much larger than the target radius, m denotes the target mass and x 

the Bjorken scaling variable ( a; = — , where is the photon virtuality ). Hence, we can 
describe DIS two step process 

1) The 7* transforms into a hadronic system well before the interaction with the target. 

2) The produced hadronic system interacts with the target. 

Gribov added two technical assumptions, which simplified the calculation 

a) The hadronic interaction is a black disc interaction. This assumption was made for a 
heavy nuclear target. In the black disc limit, the strong interaction diffractive dissociation 
channels M^p M'^p with M ^ M' can be neglected. 

b) A dispersion relation, without subtractions, can be written in terms of the variable 
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a(M,M;s) 



Figure 1: The generalized Gribov's formula for DIS. 



The resulting DIS cross section is then written as 

aem fRiM^)M^dM^ 



Here R{M'^) is defined as the ratio 



( g2 + M2 )2 
o"(e^e^ hadrons) 



(2) 



(3) 



The notation is illustrated in Fig.l where is the mass squared of the scattered hadronic 
system, V'^[M'^) = Ri^M"^) and aM'^N{s) is the cross section for the hadronic system to 
scatter off the nucleonic target. 

Assuming Eq. (d) for a^-pN and integrating Eq. (0) over we obtain a 7*p cross section 
which behaves hke ln( ^2"''''X% ) ■ We note that M^, 
one easily obtains the Gribov's bound p 



4:1711, whereas M^^^ oc s, therefore, 



■Rnc,Cln — m — 

So X 



(4) 



here R^o denotes the ratio given by Eq. (^) in the high energy limit. The logarithmic 
behaviour of the above bound is unchanged by the introduction of an arbitrary high mass 
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cutoff — < 0.1 0. A disturbing feature of Eq. is tliat it is less stringent tlian tlie 
Froissart - Martin bound, and we question wlietlier tliis is a genuine feature of DIS, or an 
artifact of our approach and the assumed input. In particular, an extra logarithmic power 
of which appears in Eq. (Q) due to the upper limit of the integration, does not appear 
explicitly in Gribov's formula. 

Attempting to clarify these problems, we present a calculation that maintains Gribov's 
basic hypothesis, which we find attractive , but gives up the black disc assumption for large 
values of M^. The physical reason why the black disc assumption cannot be correct, even 
for a very heavy nucleus, or at extremely high energies in the Froissart - Martin region, is 
simple. The quark - antiquark pair with a large mass has a small transverse size, typically 
of the order of oc -g. Being colour neutral such a pair can penetrate without interacting 
through a large target such as a heavy nucleus, or a hadron at ultra high energies. The 
attractive feature of this scenario is that at large the interaction takes place at small 
distances, or in other words, it is a hard process which can be calculated in the framework 
of perturbative QCD (pQCD). 

In section 2 we develop a general method to take into account the effect described above, 
and show that Eq. (11) should be replaced by the following relation: 



.(7 N) < j-R^C In - ln( ^ ) , (5) 



where M^(x) is the solution of the equation 

47raQ 



3R%M^{x) 



xG'^^^^^(x,M^(x)) = 1 . (6) 



xG denotes the gluon density of the target, and Rj^ the gluon correlation radius, that has 
been estimated using the HERA diffractive dissociation data in Ref. 0. Mq is a cutoff in 
mass that separates the "soft" ( lomg distance ) processes from the "hard" ( short distance ) 
ones. Clearly, the assumption that the production for < Mq is soft, whereas > Mq is 
hard, is an oversimplification. In our approach, the value of Mq is a pure phenomenological 
parameter, which we determine from a fit to the experimental DIS data at sufficiently small 
values of Q^. 

In section 3 we develop a phenomenological approach based on our general formulae of 
section 2, that allows us to match the DIS and real photoproduction data. The main idea 
underlying our approach is to parameterize the low mass region using Gribov's formula for the 
"soft" processes, while for the high contribution, the leading agln(l/x) approximation 
of pQCD is used. We show that with the choice Mq ^ 5GeV^, we are able to qualitatively 
describe the main features of the experimental data for photon nucleon interactions at high 
energy, and for most values of the photon virtuality { )■ 
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2 General formalism and a Gribov - Froissart type 
bound for DIS 



2.1 A generalization of the Gribov's formula 



As mentioned above , Gribov argued that one can use a dispersion relation with respect to 
the masses M and M' to describe the photon - hadron interaction ( see Fig.l for notation ), 
as the correlation length lc = ^ Rn, the target size. Based on this idea we can write a 
general formula for the photon - hadron interaction, 



= 37 y Q^ + M^ + M- ■ 

In the black disc approximation aiM"^, M'"^ , s) = 2 7tR%M'^ ^(M^ - M'^), which leads 
to the Gribov's formula of Eq. @. Eq. (0) enables us to separate the "soft" and "hard" 
interactions, by introducing a separation scale Mq in the integral over the masses M and M' 
in Eq. (^. Eq. can be rewritten in the form 

a{-f*N) = a'"^' + a^"'''^ , (8) 

where 

aem f^^o R{M^)M^dM^ 

3 71 J Ami ( Q + ^ ) 

Here, we have used Gribov's ideas to estimate the contribution of the "soft" processes using 
the black disc approximation. It would be preferable if the "soft" contribution could be cal- 
culated in nonperturbative QCD (npQCD ). Thus far, unfortunately, no consistent npQCD 
approach has been developed for this contribution and what we have at hand are only phe- 
nomenological parameterizations describing the "soft" hadron processes. We observe that 
experimentally -^y' decreasing rapidly with energy in nucleon - nucleon interactions [Q. 
In our context this translates to the observation that a transition from a hadronic system 
with mass M to one with mass M' is somewhat smaller than the elastic cross section which 
does not change the value of the mass. Moreover, we expect theoretically the ratio to 
decrease logarithmically in the high energy limit |@]. These observations support our idea 
that a suppression of M to M' transition can be used as the first order estimate of the "soft" 
contribution in Eq. (^. Clearly, at this stage, any description of the "soft" term in Eq. (H) 
has to be based on a model. 

An attractive feature of our approach is the introduction of the separation scale Mq. This 
allows us to use both the "soft" high energy phenomenology as well as the pQCD calculation 
for the photon - hadron interaction at high energy. 
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Accordingly, for the short distance contributions we have 



^hard 



where we can use tlie leading log(l/x) approximation of pQCD to evaluate this integral. 

2.2 The "hard" contribution to the generahzed Gribov's formula 

We wish to rewrite the formula for the "hard" DIS cross section in a form which is similar 
to Eq. (plQl). The cross section for DIS in the region of small x ( high energy ) in the leading 
log(l/x) approximation of pQCD, has the form ^ 



QCD 



where '^'^* is the wave function of the virtual photon. Although the separation between the 
"soft" and "hard" sectors is more natural in the analysis of longitudinal polarized photons 
we limit our discussion at this stage to transverse polarized photons as this gives the domi- 
nant contribution to the total cross section. The calculations pertaining to the longitudinal 
polarized photons will be published elsewhere. 

For a transverse polarized photon we have [0 

= '^Y.Zj[z' + {l-zf]Q-Kl{Qr^), (12) 

where K\ is the modified Bessel function, = Q'^z{l — z), and Nc the number of colours. 
Zf and z are the fraction of the charge, and the fraction of energy carried by the quark. r± 
denotes the transverse splitting between the quark and antiquark. crAr(x, r_L, 6^) is the cross 
section of the colour dipole of a size with the target at fixed impact parameter bf. This 
cross section is equal to 

aN{x,rj_,bt) = 2Imaei{x,rj_,bt) , (13) 

where is the elastic amplitude in the bf - representation, which is closely related to the 
scattering amplitude of the dipole at a definite value of the transfer momentum squared 
t=-q^ 

aei{x,r^,bt) = ^ I d^qte-"^'-''A{x,r^,t) . (14) 

Since high energy experimental data suggest that Rea^i <C Imaei, s-channel unitarity 
implies 

cTN{x,r^,bt) = 2{ 1 - e-^^ } , (15) 
*This subject will be further discussed in section 3. 
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with arbitrary real function f2. In the kinematic region where -C 1, for dipoles of small 



sizes r± -C Rn, this function is equal to ||10 



Q = S{h)^rlxG{x,^) . (16) 



S{bt) is the nonperturbative two gluon form factor, which normalizes as / d%tS{ht) = 1. 
From general principles of analyticity we only know its large bt behaviour, S{bt)\bt^oo 
^-2fibt^ where /i is the mass of the lightest hadron ( pion ). 

Many practical applications assume an exponential parameterization for S{bt) of the form 

1 -4 

Sih) = e ^ , (17) 

where R% is the correlation length between two gluons in the proton. For the case of 
uncorrelated gluons Rn is the hadron ( proton ) radius. 

We wish to comment on the form of Eq. (|16]), recalling the standard procedure for solving 
the DGLAP evolution equations. 

1) The first step: we introduce the moments of the parton density, 

xG{x,Q^) = — [ e"^ M(cj, Q2) rf^^ 
27rz Jc 

where the contour C is located to the right of all the singularities of the moment M{uj, Q"^). 

2) The second step: we find the solution to the DGLAP equations for the moment 

dM{uj,Q^) 



dlnQ2 

The solution is 



7(a;)M(a;,Q^) 



M{u,Q^) = M(u;, Q2) e7H inWVQg) . (19) 

Here M{u, Ql) is the nonperturbative input which is taken either from experimental data 
or from the "soft" ( model dependent ) phenomenology . 

3) The third step: we find the solution for the parton density using the inverse transform 

xGix,Q') = [ ^ e-MiA) + 7Hin{QVQg)M(^,Q2) . (20) 
Jc 27ri 



We conclude that in order to find a solution of the DGLAP equation we need to know the 
nonperturbative input M{uj, Qq) and the anomalous dimension 7(0;), which we can calculate 
in pQCD. To obtain the bt - dependence of the deep inelastic structure function we have to 
calculate the t - dependence of the imaginary part of the virtual photon Compton amplitude 
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( see Eq. (0) and Eq. (0) ). In the framework of the DGLAP evolution equations we 
have two different regions of t : (i) t < Ql and (ii) t > Q^. For the case when t > Ql, 
t defines the factorization scale and replaces Ql in Eq. (^) ( see Ref.|]n[ ). However, for 
t < the factorization scale is equal to and the only t - dependence is concentrated 
in M(a;, Qo! The factorizable form of the initial moments M(a;, Qq; t) = M{uj, Q^) F{t) is 
certainly an assumed model, but this model is reasonable for t <^ Q^. It should be stressed 
that this assumption which led to the explicit form of Eq. (^), is not needed for the large 
bt behaviour, which is the only ingredient of Eq. ( [T^ ) used for the proof of the Gribov - 
Froissart bound for DIS. 

Using Eq. ([T5|), we can distingwish between two kinematic limits that we use for our 
approximation. 

I) ^7 ^ 1 and <7N{x,r±,bt) — > Q with Q given by Eq. (p^). 

II) > 1 where a^ix, r^, bt) = 2{1 - e-^^ } ^ 2 . 



At each fixed x and r±_ the boundary between these two regions occurs at bt = bQ, which 
can be determined from the equation 

S{bo)^rlxG{x,^) = 1. (21) 
6 r_|_ 

Substituting the large bt behaviour of the form factor S{bt), one finds 

60 = ;^ ln[rlxG(x, -^)] . (22) 

In the kinematic region II we rewrite Eq. ( pi]) in the form of Eq. (0). This is a very 
simple task once we recall that 

Q^TT' ^o(fcr^) , (23) 

or 

kldki f kldk2 

Using the simple form of aN{x,rperp, bt) = 2 in region II, one can integrate Eq. (|ll|) over k2 



■,QK,{Qr.)r = J ^^Mhr.) J Mhr.) ■ (24) 



9 k 

and z introducing a new variable M = \ - We obtain 



2(1-2)- 



2lZf / , T7^r—T7^2^ i dhf (25) 



At first sight it appears that there is no natural cutoff in the above integration, but as 
we shall see later this is not so. 
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2.3 The unitarity bound on the photon cross section 



To obtain the unitarity bound for the total cross section of the photon - nucleon interaction, 
we use the decomposition of Eq. (|^) and Gribov's estimates for 0"*°-^* given in Eq. (P). For 
c"M2iv(s) in Eq. (|^) we can apply the Froissart - Martin bound of Eq. (0), since it is a 
typical hadronic ( on mass shell ) cross section. Note, that the Froissart - Martin bound is 
a high energy limit for which the Gribov black disc assumption is perfectly adequate since 
the diffractive ( M ^ M' ) channels are suppressed relative to the elastic ( M = M' ) 
channel. We can evaluate cr^"^'^ in Eq. using the inequality a^"-^'^ < a^'"^ where a^'"^ 
is determined by Eq. p5|). We wish to explore further the nature of the upper limit of the 
integration with respect to in Eq. (|25|) . To do this we return to Eq. (^Tj). This equation 
has no solution if 



S{bt = 0) 



r_|_ xG{x, —^) < 1 . 



(26) 



Therefore, the main contribution to cr'^^^ in Eq. ([TTD comes from the kinematic region II 

is a solution of the equation 



with > Tnfx 



where r^i^x) 



S{bt 

Using the notation S{bt = 0) 
familiar form 



0) — — rl{x) xG{x, 



(27) 



TTRi 



see Eq. (0)), Eq. ( p7|) can be rewritten in a more 



— r2-ro(x) xG(x, — 



3R 



N 



Due to the uncertainty principle 



oc 



r^{x)' 



1 



1 . 



As the integration]] over z in Eq. (^TJ) is convergent in the limit (JN{x,rpf,rp,bt) 
Eq. ( p4D and Eq. (|^) ), we can safely put z = | in Eq. (|^)|, and rewrite Eq. 
equation for the upper limit of the integration over 



(28) 



(29) 

2 ( see 
81) as an 



Anas 



xG{x,M'^{x)) 



3R% M^{x) 
Collecting all contributions we find 



a{j*N) < ^ { C In^ - / 

3 TT So Jiml 



+ 2R^bl ln( 



Q2 + M'^{x) 



Mi 



(30) 



)}■ (31) 



This equation provides an improved Gribov - Froissart bound for the photon - hadron total 
cross section, in place of the less restrictive one given in Eq. (^). For very small x Eq. ( pO]) 

^We will comment on the z 
tJt means that r — 
formula of Eq. (|^) . 



- integration in the "hard" cross section for ^ 1 later, 
oc This fact justifies our main input of a separation scale ( Mo ) in Gribov's 
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leads to M'^(x) h? exp{^ a\xi{l / x)) ^ where A is the QCD scale and the constant a has 
been calculated in Ref . ||1 If . This implies that for very small x, the bound for (7(7* A^) at 
g2 < M2(x) is 

o{l*N) < ^ 2 [ JaHl/x) - In(gVA') ] C In^ - . (32) 

In the ultra small x limit we obtain that cr(7*A^) < C ( ln(- )t where C contains all relevant 
constants. 



2.4 Numerical estimates for the behaviour of cr(7*A^) 



We can use Eq. ([HT| ) also to make numerical estimates of the high energy behaviour of the 
DIS total cross section. For this purpose we rewrite Eq. (|3l|) and attempt to evaluate au'^ n 
of Eq. (P) rather than using its high energy bound. To this end we assume 



a. 



37r 



r soft 



hadron 



4m,r 



+ M2)2 



+ MS 



(33) 



where crl°Jl.^^{s) is a typical cross section for a meson - nucleon interaction. To obtain an 
estimate, we take crladroni^) = k ( '^(^'''P) + cr(7r~p) ), and use the Donnachie - Landshoff 



parameterization [14| for its energy behaviour. 





Figure 2: Solutions of Eq. (|2T]) {bl{x) ) and of Eq. (jS^) (M'^ix)), using the GRV parame- 
terization for the gluon structure function. 



Solutions of Eq. (|TD ( i.e. feg(x) ) and of Eq. (ID ( i.e. M'^{x) ) are plotted in Fig.2, 
using the GRV parameterization ||T^ for the gluon density. In solving Eq. ([21|), we have used 
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Eq. (|r^), with R% = lOGeV^. The value of R% is derived from HERA data on diffractive 
production of vector mesons in DIS, and from the high energy phenomenology for "soft" 
processes ( see Ref.^ for details ). We need to calculate 6q at fixed x for all values of < 
M^(x). However, since the integral over is logarithmic, we can evaluate its contribution 
at an average M^. is determined from the relation J^2 qt^j^ = 1^2^^^ q^+1p ' "which 



gives = y'(Q2 _^ M'^{x) + M^)-Q'^. hi at this mass value is plotted in Fig.2. We 
consider the values, given in Fig.2 , to be more relevant at presently accessible energies, than 
the highly asymptotic Froissart - like estimates. 

Fig. 3 shows the energy dependence of the r.h.s. of Eq. (^) together with the experimental 
data. The values of M'^{x) ( see Fig. 2b ) show that we can trust our estimate, given by 



Eq. (33), only for x < 10~^. This means that we can compare our bound only with the 
available experimental data at relatively small values of Q^. We plot the data and our 
estimates from Eq. (|33D only at ~ 0, since for quasi real photoproduction we reach 
the smallest values of x. One should note that Eq. (|33|) was proven only at very small x, 
where we can neglect the contribution from the kinematic region where f2 <^ 1 ( region I ). 
Actually, in the HERA kinematic region at all available and W, we cannot neglect the 
contribution from region I. On the other hand, Fig.3 shows that our bound is rather close 
to the experimental data for real photoproduction. This suggests that Eq. (^) reflects the 
main physics for the HERA kinematic region. We expect that the inclusion of the "hard" 
contribution from kinematic region I, will improve our estimates, but will not produce a 
dramatic change. 



3 Matching of the "soft" and "hard" processes in DIS 
3.1 General description 

In the following we develop a phenomenological approach to describe DIS at all values of 
based on the separation of the "soft" and "hard" interactions, in the framework of the Gribov 
formula ( see Eq. (0) and Eq. (H) ). This approach provides a relatively simple description, 
in which the mass integration with < Mq is controlled by the "soft" interaction. For 
> Mq we are dealing with a "hard" interaction, which we can treat in pQCD. For a^°-^^ 
in Eq. (^) we assume that the "soft" high energy strong interaction suppresses the diffractive 
dissociation of a hadron state of mass a M to a hadron state with a different mass M'. This 
property is true, for example, in the additive constituent quark model, where the interaction 
of the hadron can be reduced to an interaction of the quarks, namely, only the first diagram 
of Fig. 4 ( with M = M' ), contributes. In a different context this is also a consequence of 



the implementation of screening corrections |Tl|. With this assumption, we can rewrite the 
"soft" contribution in the form 
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Figure 3: The energy dependence of the r.h.s. of Eq. (|33D together with the experimental 
data. 
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p 



p 



p 



p 



Figure 4: The diagrams for two gluon exchange modeL 



For cr^'^^'^ we use the general formula of Eq. (|ToD, where (t'^*"^(M^, M'^, s) is rewritten in 
terms of the gluon - nucleon interaction in the framework of a two gluon exchange model, 
shown in Fig.4. This model is certainly correct in the region of very small x {as ln(l/a;) > 1) 
and large (asln(M^/A^) > 1), and it also reflects the main properties of the QCD 
interaction outside this particular kinematic region. 

We shall discuss the assumptions made for cr'^"-^* and Eq. (|3^ later. Prior to that we 
wish to specify the well known "hard" contribution. This will also be instructive for our 
discussion of the "soft" contribution in section 3.3. 

3.2 The two gluon exchange model 

Eq. ( ]TTD is the basic formula for the two gluon exchange model, where we use the following 
representation for a{x,r±) = J (Pbt(7N{x,r±,bt) 




(35) 



One can easily see that the two terms in Eq. (^) just reflect the two diagrams in Fig.4. 



Substituting Eq. ( ^5]) in Eq. (|ri|) and using Eq. (^4]) , we obtain 




(36) 
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In Eq. (|36| ) we integrate over the angle between ki± and l± and introduce a new variable M 

h.2 k2 72 

~ zil-z)' ~ zil-z)' ~ zil-z)- ^^^^ 

The physical meaning of M is clear. Indeed, on the average 

\M' - M"\ = - ih^-^J±l\ = • + =<M'>. 

' ' z{l-z) z{l-z) ^ z{l-z) 

In terms of the new variables, Eq. ( p6D has the form 

^-^^[^A^^i^^m JKQ"-^ (38) 

^ - ^ Q' + M' - ^ ^^^^^ 

+ ^(Q2 + M2 + M2)2 - 4M2M2 

Eq. d 

iVc Ef 



Note that Eq. (|38|) applies in the region M2 > Mo2 and that R{M'^) replaces i?c 



Since z{l — z) = Eq. (|38|) can be rewritten in the form 



4772 74m. M4 Jui + JqI ' M^' L I^i 

, M2 - g2 Q2 ^ j\^2 _ j^^2 

{ A^2 ^ n2 + / - } ■ (39) 

+ ^ V(Q2 + M2 + M2)2 - 4M2M2 

Recalling that a (ZD = as(/i)^, where as(Q^) ^^^(x, Q^) = as(/2) 0(/2) ^/2 ^ 

M2 



Ref. |Tl| for details), we obtain, in the limit -r^j <^ 1, that 



37r Af2 g2 _^ ^2 y^Qg 4. ' ^ ' 4 ^ ^ ^ 

^ M2 - g2 ^ g2 ^ j^^2 _ j^^2 ^ 

M2 + g2 y^(Q2 + ]Vf2 + M2)2 _ 4M2M2 

This equation is our master formula for the evaluation of the "hard" contribution ( with 
light quarks ) to the total photon - nucleon cross section. Qq is the starting value of the 
gluon virtuality for the DGLAP evolutuion equations. 

We wish to add two comments concerning the above master equation. 
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1) z has the same value before and after the interaction. This is a direct manifestation 
of the leading log (1/x) approximation in which we only take into account contributions of 
the form ( ag ln(l/x) )". To understand this we compare the time of the interaction of the 
qq - pair with the target ( Tj ) to the life time of the virtual photon fluctuating into a qq - 
pair ( r^* ). According to the uncertainty principle 

r,, ^ J- = I I I = - (41) 

where ki and k2 are the four momenta of quark (antiquark) (see Fig.4). An estimate of 
the interaction time can be obtained from the typical time for the emission of a gluon with 
momentum / from the quark ki, ( i.e., see the second diagram in Fig.4 ). 

~ Ifci^ - k\_ - iJ ^ 'Mx _ iL _ ' ^^^^ 



a 



where « = ^ and z' = z — a. In the leading log(l/x) approximation we have a <^ z and 
hence 

r, ^ ^ « T,. . (43) 
'± 

Eq. ( |43| ) shows explicitly that the processes of a virtual and/or real photon scattering off a 
target can be described as a two stage process: initially the photon decays into a. qq - pair 
and only a long time after that the pair interacts with the target. This is why such a process 
can be described by the Gribov's formula. 

2) An important feature of Eq. (^) is the fact that we integrate over z using z{l — z) = 

■j^ . The integral over is logarithmic and the typical values of are of the order 

of ~ M'^exp{—^) where 7 is the anomalous dimension. Since we are interested in 
the region of small x where 7 is rather large, one can see that the typical values of z, 
which are essential in our integration, turn out to be of the order of unity. As a result the 



contribution of the aligned jet model ( AJM ) [^] is small at least at high energies. Actually, 
an explanation showing that the emission of many gluons suppresses the nonperturbative 
AJM - type configurations has been given by Gribov and Lipatov more than two decade 



ago. They showed that the emission of multi gluons generates a Sudakov form factor which 
suppresses the AJM - like contributions. Note that the result of the z integration is very 
important for the understanding of our approach since it justifies our idea that is a good 
measure of the typical distances involved in the process. Indeed, oc holds only for 
^ ~ 1 ( see Eq. (|2^) ). These estimates have been made for the total cross section. For 
exclusive channels the situation is quite different. For example, for the diffractive dissociation 
in DIS, induced by a transverse polarized photon, the typical distances are rather large and 
do not depend on the value of the mass in the Gribov's formula. 

For heavy quarks, the diagrams of Fig.4 give 
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Kara «e„^. 2 ^ ^^^(M^) dM^ M\ 

{ — + — z — I 

where i?*3<9 jg ^j^g heavy quark contribution to the ratio in Eq. (|^), and mq is the mass 
of the heavy quark. In Eq. (^) we have assumed that the quark is heavy enough to be 
described in pQCD without any contribution of the "soft" processes. In the above formulae 
X = , and W is the energy of the photon - nucleon interaction. 



3.3 A model for the "soft" interactions 



Our model for the "soft" interactions is based on Eq. (|3^) . We observe that our discussion 
on the time structure of the photon - hadron interaction does not depend on any specific 
properties of QCD, and can be considered, therefore, as a main feature of the parton model 
approach to high energy photon induced interactions ( see Ref. ). It means that, the 
approximation z = z' vn Fig.l applies also for the "soft" interactions of our model. 



In the parton model [13|, as well as in the high energy phenomenology |Ti|, the "soft" 



processes at high energy can be described by the exchange of a soft Pomeron which has the 



property of Regge factorization |15 



where ap{t) is the Pomeron trajectory and gp{M,M') and Gp are vertices of the Pomeron 
interaction with the quark - antiquark pair and with the proton respectively. Substituting 
this equation in Eq. ([TT| ) one can see that only the photon wave function depends on z. This 
integral is convergent and z is typically about unity also for a photon with large virtuality. 
Coming back to Eq. (|37D one can see that M' could be different from M only if k2± ^ ki± 
( or k2± <^ ki± ). On the other hand, in the parton model the typical transverse momenta 
( l± in Fig. 4 ) are of the order of the "soft" scale, i.e. about 1 GeV. Therefore, in a photon 
- hadron interaction we expect that the value of M' cannot be much larger than the value 
of M. The success of the additive quark model ( AQM ) ( see the first diagram in Fig. 4 ) 
in the description of the high energy scattering ( see Ref. |14[ ) supports our assumption 
that the difference ( M' — M ) is much smaller than the transverse momentum scale for the 
Pomeron. 

We can, therefore, rewrite the general Gribov formula of Eq. (0) in the form 
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aem f^^o R{M'^)M'^dM^ . , f^^ MdAM gp{M,M') + 

-(yM^NKs) 



37r 74m. (Q2 + M2)2 ™ ^ J^M M2 gp{M,M) {M + AMf + Q2 

= ^ (Q2 + M2)2^*^^^(^) ' (45) 

where AM = M' - M. 

Note that the AM integration is particular to photon induced reactions and is, obviously 
absent in the case of hadron - hadron collisions. This correction, as well as other relevant 
corrections, are absorbed in our definition of aM^n in Eq- (HD- 

In light of the above discussion, we choose Eq. (Q) and /or Eq. ( ^Sj) as the master formula 
for our description of the "soft" contribution. 



The quantities appearing in Eq. ( p4D which need to be specified are the ratio R{M^), and 
the interaction cross section of a hadronic system with mass M with the target ( (Tm^n{s) in 
Eq. (|3^ ) ). Although, there is experimental data for R^M"^) ( see Ref. |]T^ ), we have used the 
parameterization of Ref. ||T8[| for R[M'^) |. The two main ingredients of this parameterization. 



reproducing the experimental data, are the resonance contribution and the background, 
which at high masses approaches the QCD result ( R^M"^) — > R^o )■ We approximate ^m^n 
in the resonance region by 

aR = fi;i(o-(7r+p) + a{n^+p)) (46) 



where we use the Donnachie - Landshoff Reggeon parameterization [0 for the cross section 
of pion - proton interaction, 

aji = k[A (A)"p(o)-i + B (A)"«(o)-i] ^ (47) 

So So 

with a Pomeron and Reggeon trajectory intercepts of ap(0) = 1.079 and a_R,(0) = 0.55 ; sq = 1 
GeV"^; A = 13.7 mb and B = 31.9 mb. We introduce a rescaling constant k as a parameter in 
Eq. ( PB[ ) and Eq. (0). k = 1 corresponds to a simplified AQM where the integrant is taken 
at AM = and first order corrections of and where mo is the typical scale of the 
soft Pomeron, are neglected. This is discussed further in section 3.4. 

To describe the interaction of the background contribution in R{M'^) , we need to deter- 
mine the correct energy variable for the interaction of the hadronic system of mass M with 
the target. We suggest to replace the variable ^ by where xm = In other words. 



we replace Eq. ( ^T]) by 



ft; [A' (^)ap(0)-l ^ B' (^)"«(0)-l] , (4J 



Xm Xm 



^We thank E. Gurvich for drawing our attention to Ref. 1 18 
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where A' = 13.1 mb and B' = 41.08 mb. The values of A' and B' were chosen so that 
Eq. (^) for the p meson - proton interaction is vahd. The variable xm is not unique, but 
it correctly describes the high energy interaction, in all parton - like models in the region of 
large mass M ( see Ref . ITlI] for example ) . 



3.4 Comparison with the experimental data 

To compare with the experimental data, we calculate the total cross section of the photon - 
proton interaction using the following formula 

cr(7*p) = a^"f\Eq. (34), Eq. (H) ) + a^^^^Eq. (0)) + ^-'^(i^g. (H ) . (49) 

Eq. (^) depends on two free parameters: k and Mq and the input gluon distribution 
xG{x, ^). In Fig. 5 we plot cr{'y*p) as a function of W'^ for various values of and compare 
with the relevant experimental data. In our calculations we have used the GRV parameteri- 
zation |T^ in higher order of pQCD (GRVHO). We have chosen the parameters k = 0.6 and 
Mq = 5 GeV^ so as to obtain a reasonable reproduction of the data. 

As can be seen from Fig. 5, we obtain a good fit fot the low data over the entire W 
range. However, at higher than a few GeV^, our description of the data is deficient in as 
much as we are not in agreement with the low energy experimental points and our predicted 
high energy behaviour is steeper than the data. Below we elaborate on these features as well 
as on some important details of the suggested model. 

1. With only two free parameters, we manage to reproduce the energy dependence for 
real photoproduction and DIS with < 8GeV^ cross sections. This is shown in Fig. 6 
where we compare the recent high energy - low ZEUS data pS) with the predictions of 



our model. Our results compare favorably with the Donnachie-L andshoff [P^ and GRV 
parameterizations. 

2. The fact that we are unable to reproduce the low energy behaviour of the higher 

( X > 10"^ ) data is not surprizing. These data correspond to higher x values for which 
the two gluon approach to DIS is insufficient. In Fig. 5 we show the line corresponding x = 
0.01 which illustrates this point. Clearly for high x ( low W ) one should also include the 
contribution coming from the evolution of the quark distribution as a part of the pQCD 
description of DIS. We discuss below how to expand our formalism so as to include this input 
as well. The fit can also be improved by the utilization of other input parton distributions 
which are less steep than GRVHO in the small x limit. 

3. An unexpected feature of our results is that we require a value of k = 0.6 to rescale 
the AQM estimate of the "soft" contribution. As we have noted, this value of k, refiects 
the need to integrate over AM which is particular to photoproduction and DIS and does 
not appear in hadron - hadron scattering, where the AQM has been checked experimentally. 
Our result is different from VDM |T^, where in order to describe the experimental data, 
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Figure 5: The comparison of the experiment data ( see Ref.[p^ and references therein) for 
o"(7*p) and our calculation, using Eq. (^^. The diagonal line indicates the boundary for 
X = 10-2. 
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one has to assume that the vector meson - nucleon cross section is bigger than the AQM 
estimates. This difference arises from the background contribution that is neglected in the 
VDM approach. It should be stressed that such a contribution, which is included in our 
parameterization of -R(M^), is needed to reproduce the - dependence, which is much 
smoother than the VDM prediction. In addition to uncertainties from the AM integration 
we also see at least two reasons leading to a value of k smaller than unity. Our evaluation 
suggests that using the approximate formula of Ref. |TB| we overestimate the experimental 
data by about 10% ( k 0.9 from this source ). The second source is the need for shadowing 
corrections (SC ). SC definitely diminish the value of the cross section. We can estimate 
the SC from the value of the diffractive dissociation cross section, using the AGK cutting 
rules [^, which relates the SC to the total cross section, namely, Acr'^'" = o"^^, where 
a = a^'^^^ — /S.a^^ . The experimental value of the diffractive cross section is about 15%, 
both from real photoproduction and from pion - proton interaction. These two sources 
suggest a value of k ~ 0.7 to which an additional 15% reduction is added due to the AM 
integration. 



4. The inclusion of heavy quark contribution is important. In Fig. 7 we plot the ratio 
a(charm quarks) ^ functiou of energy at different values of Q^. This ratio is 

u(hght quarks) ^ 

rather small for real photoproduction, reaching 5% at high energies. For large values of 
the ratio approaches 30 - 35% at high energies W > ?>Q GeV. 
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Figure 7: The ratio = quZS function of W for different values of 
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ratio Rg 
toi?f ^ 



In general, our model provides a very simple approach to incorporate the "soft" and 
"hard" components of photoproduction and DIS. In the detailed fit of the data we observe 
two features that we consider to be rather general. 

1. The high energy "hard" contribution turns out to be sizable, even for = 0. To 
illustrate, how important the "hard" contribution is in our formalism, we plot in Fig. 8 the 

= lTo7tl^Z) - ^^^^ = 0' ^5 ~ 1% at W = 10 GeV and it grows 

15% at W = 300 GeV ( s = W'^ ). This increase is sufficiently large, that it alone 
may account for the experimentally observed increase in the energy dependence of the total 
cross section, for real photoproduction. In other words, it is possible to fit the experimental 
data assuming that the "soft" Pomeron ( see Eq. (^) ) has an intercept which is equal to 
unity ( ap{0) = 1 ). 

Our result suggests a possible and interesting interpretation for the origin of the exper- 
imentally observed increase of the total cross sections for hadron - hadron interactions. As 
for each hadron we also have a contribution of the "hard" process to the total cross section, 
due to the probability that two quarks approach each other at a sufficiently small distance. 
The probability for this "hard" process is controlled by the respective wave functions of the 
interacting hadrons. 
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2. We find contamination of the "hard" processes by the "soft" ones. For example, at 



soft 



(7*P) 



= 15GeV^ which is usually considered a large value for Q^, the ratio Rfj 
changes from 1 at W = 30 GeV to 0.2 at W = 300 GeV. Even at 
is about 0.12 at W = 100 GeV. The above results, lead one to view DIS in a new hght 



(t'""''*(7*p) 

65GeV^, Rfj 
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and provide a basis for a better understanding of what is meant by small distances or high 
photon virtualities. 



3.5 A generalization of our model to larger x 

The generalization is based on Eq. pO|), utilizing the fact that for the DGLAP |^ evolution 



equation the following equation holds in the region of small x ( see for example ) 

QpDGLAP 



2 ^"S ^r^DGLAP 



91n(QVA2) 97r 



^^^^^^^(x,^^) . (50) 



QpDGLAP 



We suggest replacing asxG{x, Q ) in Eq. (^D|) by gin(Q2/A2) ^"^^ ^^^^ generalized formula 
for DIS, even in the region of x not very small. Note, that for x not too small, we obtain 



see Ref. for example) a more general formula for Qi^iQi/i^i\ than Eq. (|50|) which includes 



ain(QVA^) 

the quark densities. After doing so, Eq. (HO) reduces to the form 



Hard , 2 T ^(^') T dF^<^^^P{x, ^) 



M2 + Q2 J(g2 + M2 + M2)2 _ 4M2M2 



Integrating Eq. (^) by parts one obtains 

roc i?(M2) roo m2 

^ "^"^ A/,? Q2 + M2 70 M4 "^2 4 J l^^J 



^hard 



q Q2 + j\^2 7o ^4 

M2 - Q2 Q2 ^ jg2 _ j^^2 4Q2M2M2 



+ 7(g2 + M2 + M2)2 _ 4M2M2 [ y(g2 + M2 + M2)2 _ 4 M2 M2 



{ 

Actually, a formula of the same type as Eq. (^) was first suggested by Badelek and Kwiecin- 
ski but using our formalism we obtain quite a different result. We can consider Eq. (p2|) 
as a generalization of the Badelek - Kwiecinski approach. In addition to the resonances we 
have also included the background contribution, and obtain the contributions of both "soft" 
and "hard" processes by integrating over M"^ and M"^ in Eq. (|5^). 



Numerical results pertaining to Eq. (B2]) will be published separately. 



4 Conclusions 



We have achieved two goals in this paper: 
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1. We provide an explanation of how and why the short distance ( "hard" ) interaction, 
calculable in pQCD, provides a mass cutoff in the Gribov's formula for photon - hadron col- 
lisions. We have shown that the Gribov bound ( see Eq. ) given in Ref.0), overestimates 
the photon - hadron cross section, and should be replaced by a more restrictive limit , as 
derived in this paper ( see Eq. (PT|),Eq. and Eq. ). At fixed as IV — >• cxd our 
bound is cr(7*p) < C'iln-)'^. 



2. We developed a simple model which consists of two contributions: "soft" and "hard". 
The "soft" term describes the long distance contribution, while the "hard" term is related 
only to the short distance interaction controlled by pQCD ( and the DGLAP evolution 



equation |^ ). This simple model with only two parameters provides a good description 
of the available experimental data over a wide range of W and < SGeV"^. We have 
suggested a technique of how to improve the high Q"^ results at sufficiently small values of 
energy W. 

Examining our model we found two interesting features that may be more general: 

a) Short distance effects contribute even at = for high energies. The contribution 
is sufficient to explain the energy rise of the total cross section for photoproduction, which 
has been interpreted as an argument that the "soft" Pomeron has an intercept larger than 1 
( see Ref . ||14|| ). This result encourages us to reconsider this widely held explanation, and to 



estimate the contributions of the "hard" processes to the growth of hadron - hadron cross 
sections, with increasing energy. 

b) The long distance processes contribute to the total cross section even at rather large 
values of Q^. For example, at = 65GeV^ and W = 100 GeV^ they are responsible for 
10% of the total cross section. This observation can be very important for understanding 
the energy dependence, as well as the value of the cross sections of other processes such as 
diffractive dissociation, inclusive production etc. We propose to examine these processes in 
the near future using the same approach. 

Our approach is not in contradiction with the usual description of "hard" processes, 
based on the DGLAP evolution equations with initial nonperturbative parton densities at 

= Qo- However, we significally enlarged the region of applicability of such an approach, 
noting that the quark - antiquark pair with large mass can be treated in pQCD even at 

= 0. It allowed us to separate the nonperturbative contribution in a different way than 
usually done and to calculate a part of the initial parton densities at = Ql in pQCD. 

In general, the model suggested allows one to discuss the interface between long and 
short distance processes, not only on the qualitative level but also on a quantitative one. 
We are of the opinion that our model incorporates what is known, both theoretically and 
phenomenologically about "soft" and "hard" physics, and provides a method to estimate the 
different contributions to a variety of processes. It also allows one to specify the kinematic 
region where the "hard" contribution dominates, and to calculate it within the framework 
of pQCD. 
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